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Abstract 

We present simple explicit estimates for the apparent reaction rate constant 
for threemolecular reactions. For small concentrations and d > 1, it depends 
only on the diffusion coefficients and sizes of the particles. For small concen- 
trations and d < 1, it is also time - dependent. For large concentrations, it 
gains the dependence on concentrations. 
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I. INTRODUCTION 



Threemolecular processes are important for catalysis ||T| , where two particles, say A and 
B, react in presence of the third substance, "catalytic cite" — C. Possible applications to 
coagulation were also concerned 0]. 

To our knowledge, the first attempt to obtain a simple theoretical description of the 
kinetics of the threemolecular reaction^ 

A + B + C^ Product, (1) 

where kt is the "chemical" reaction rate constant, which determines the reaction rate per 
one threeparticle encounter, was made in @J^, where the heuristic mean-field like analysis 
of the reaction kinetics for large C — sites was presented. Recently this subject gained more 
interest ^M- 

The theory of the himolecular diffusion-controlled reactions 

A + B^ Product (2) 

started with the mean-field type approach in ci = 3 P|. The most important result of this 
theory is the fundamental Smoluchowsky reaction rate constant, ksmoi, which determines 
the apparent rate constant kapp and, therefore, the mean reaction rate 

= ~kappC aC B} (3) 

where Ca and Cb are the mean concentrations. Frequently kappoheys the "inverse resistance 

law" iig 

111 /.N 

+ Z. ■ (4) 



kapp kfransp k chem 

For the reaction (0) in 3(i the "transport" constant ktransp = ksmoi and the "chemical" 
constant kchem = kb- First attempts to improve the Smoluchowsky theory concerned correc- 
tions to ksmoi and it took a long time, starting from the early works |TTHl^, to recognize 



The Product can contain C and/or A or B. 
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that fluctuation effects, wfiicfi determine tlie long-time asymptotic for the concentration de- 
pendencies, can be important. Contrary to it, in threemolecular reactions theory, the most 
of the recent works are concentrated on the long time behavior (See, however, [^] where the 
whole time domain was studied by means of a new elegant method) and/or = 1 in spite of 
the fact that complete mean-field like theory is not finished and the analog of ksmoi "was not 
calculated. Particularly, @J^,0,§ stress that for the long-time kinetics of the particular type 
of the reaction (|l]), with A = B = C , i.e. for A + A + A I A with / < 3, the fluctuation 
effects are not decisive and that the conce ntration of A, Ca, for d > 1 is governed by 



Ca oc 1/Vt, (5) 
and for d = 1 contains "logarithmic corrections" 



CA^^/{lnt)/t. (6) 

In our paper we present simple estimates for the Smoluchowsky - like rate constant, kf^, 
which determines for threemolecular reactions the apparent rate constant in the equation 

dCA{t) 



dt 



-kaj,p6ACA{t)CB{t)Cc{t), (7) 



where 6a = the number of A particles, which enter the reaction, minus the number of A 
particles which the Product contains, for all d. We show that: {i) For small C — sites 
Eq. (0) leads to the dependence Eq. for (i > 1. For d = 1 we regain the result, Eq (^, 
since kth becomes time dependent. We stress that in this case the result is still mean field 



one, since the fluctuations are not incorporated (see [T^ for more detailed analysis of the 
fluctuation effects on the k - particle reactions with different types of initial distributions). 
However, even for d > 1, kapp is renormalized by the transport processes and, therefore, it is 
not equal to kf. (li) For large C — sites kapp in the Eq. (|^ becomes a function of the mean 
concentrations and, threfore, Eqs. (|^,§) are not valid. 

In order to understand the main scaling laws for kth we consider first the simplified hoping 
picture with correlations in the reaction zone and without correlations in larger distance. For 



important limiting cases we present also the more traditional, diffusion - reaction equation 
approach. To make clear the ideas we start with the bimolecular reaction (^. 



II. BIMOLECULAR REACTIONS 

A. Smoluchowsky theory 

According to the Smoluchowsky theory, the concentration of B particles in the distance 
r from the center of the A particle, CB{r,t), is governed by the diffusion equation 

= DABV^Cn{r,t), (8) 

where Dab = Da + Db, with the boundary conditions 

$ = SdDAB^CBir, t) \r=R^^ = hCBiRAB, t) , (9) 

CB{r = oo,t) = Cb, (10) 

where Rab = Ra + Rb and $ is a flux of B particles through the d — dimensional sphere, 
|r| = Rab, with the surface Sd- The reaction rate is equal to the quasi steady state value of 
$ multiplied by Ca- For c? = 3 it leads to the Eq. with kchem = h and ktransp = ksmou 

At^RDab for d = 3 



^Smol 



for = 2 (11) 



In 



DAB/T^t for d = 1. 

For t ^ R\q/D in d < 3 and for ^ ksmoi for d >3, kapp = ksmoi- In = 3, the deviation 
of Csir, t) from the limiting value, Cb, decreases proportional to Rab/I^I and, therefore, the 
size, L, of the "reaction zone" , where the distributions of A and B are correlated, is of the 
order of the reaction radius Rab- Contrary to d = 3, in d = 1 the correlated region grows 



proportionally to y/DABt- The d = 2 case is marginal. For smaller dimensions, d < 2, the 



diffusion is recurrent, the space exploration is compact |18], which means that the number 



of the returns of the diffusing particle to the origin tends to infinity if t oo. Therofore, 



the volume of the reaction zone equals to [IS 



(12) 



where 



L oc Rab for (i > 2 and L ex. J DAst for < 2 



(13) 



B. Hoping model 

Let us now consider the simplified picture of the reaction (H), where the correlations for 
A and B are present only in the reaction zone and the exchange of A{B) particles between 
the reaction zone and the non disturbed region is a "one step" process with the frequency, 
T^A{B) OC Da{b)/L'^, which equals to the minimum eigenvalue for the corresponding diffusion 
problem. In this section, we also assume that the concentration of the reacting particles is 
small, so we can consider only pairs AB and neglect the configurations with more particles 
B{A) in the reaction zone of A{B). Therefore, the reaction rate is proportional to the 
number of AB pairs. Nab, and the balance equation reads 

^ = UAB {QVCaCb - Nab) - h^, (14) 

where V is the total volume of the system and uab = i^a + i^b- The first term in the r.h.s. 
of Eq. (14) is the rate of jumps into the reaction zone, the second term is is the rate of 
jumps from the reaction zone and the third term is the reaction rate. The quasi steady state 
solution of Eq. ( p!4D for the reaction rate leads to the ressult pesented by Eq. (|^) and Eq. (^ 
with kchem = kf, and hransp = i^AB^- Taking advantage of the Eqs. (pTD-(p!^), we conclude 
that 

l^AB^ = ksmol (15) 

with the accuracy of the insignificant numerical multipliers for d = 1 and (i > 3, and with 
the accuracy of oc log(t) corrections for d = 2. 



III. THREEMOLECULAR REACTIONS. 



A. Probability distributions. 

Consider the reaction (|l]) with Dq = 0. In order to take thremolecular correlations into 
account we write down the master equation for the joint probability, P{nA,nB), to obtain 
Ha particles A and particles B in the reaction zone, near a C — site, 

dP{nA,nB) 
dt 

The diffusion and reaction terms are equal to 



I react • 

(16) 



hiff = - {i^AnA + i^BnB + i'a^Ca + ^b^Cb) P{nA, ub) + 

{ua + l)z/A-P(^A + 1, nB) + {ub + l)uBP{nA, riB + 1) 

+ VA^CAP{nA - 1, Ub) + UBnCBP{nA, UB-l), (17) 

Ireaet = - R{nA, Ub) P {ua, Ub) + 

R{nA + l,nB + l)P{nA + l,nB + l), (18) 

where R{nA, ub) is the reaction rate for the reaction zone with the given numbers of A and 
B. We assume that the local reaction rate is equal to the product of A and B concentrations 
in the reaction zone, 

R{nA,nB) = kt ^^2^ ■ (19) 



B. Hoping model 

1. Small A and B concentrations (small C — sites) 

For small A and B concentrations, Ca{Cb)^ « 1, we consider only triples, the particles 
C, which have A and B in the reaction zone, and pairs, the particles C which have A or 
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B in the reaction zone. The mean number of C particles is Nc = CcV; of CA pares — 
NcA = NcP{l, 0}; and of CB pares — Ncb = NcP{0, 1}. The reaction rate is proportional 
to the number of the triples, Nqab = NcP{l, 1}. From the Eqs. ([1^ )-(p!9D we obtain 

= l^Ai^VCcCA{B) — NcA(B)) — 

pb{^Nca{b)Cb{A) + Ncab). (20) 
= M^VCcCb - Ncb) + 

va{^VCcCa - Nca) - hn-^NcBA. (21) 

When ua » i^b^Cb and ub » ^a^Ca, the steady state state solution of the 
Eqs. (13) '(0) Isads to the Eq. (^ for the reaction rate and to Eq. (P for kapp with kchem = h 
and ktransp = i^AB^'^- Taking advantage of the Eqs. (|TT|)- ([r3|) , we conclude that with the 
accuracy of the logarithmic corrections for d = 1 the threemolecular analog of the Smolu- 
chowsky constant equals to0 

kth oc VAB^'' « R^'^-^Dab. (22) 



When one of the diffusion coefficients, say Db, tends to zero, Eqs. (^Dj),(^T]) lead to 

= -vb^SaCcCb oc -ksmol^ACcCB- (23) 



2. Large A or B concentrations (large C — sites) 

For Ca{b)^ » 1, the analysis of the Eqs. (0)-(0) for the averaged A{B) concentrations 
in the reaction zone, C'xIb) = E5^=oE^=o -P('^a, nB)?2A{B)/f^, leads to 



%a = - 4',',,) - t.^^. (24) 



^ Note, however, that in this case the reaction on the catalytic surface with the typical time 

/ 



— 1 (2) 

'^surf ^ Ra.bCa b(^A + Db) can become the limiting stage of the reaction 



Note that the reaction term (the last term in the right hand side of the Eq. (|^)) in this 
regime decouples into the product of the concentrations. For the quasi steady state Eq. (^ 
) reduces to an algebraic second order equation for the reaction rate. When ktCs << i^a^, 
it leads to the Eq. (|^) with 

When ktCs » i^a^ and kfCA » i^b^, the steady state Eq. (p^ leads to 

dCA 



dt 



QCcmm{uACA,iyBCB). (26) 



C. Diffusion appoach 

1. Small A and B concentrations (small C — sites) 

For the small concentration limit we propose the following extension of the Smolu- 
chowsky approach to threemolecular reactions. We place the origin into the cen- 
ter of an immovable particle C and determine the the conditional density of triples, 
CAB{rA,rB,t)d'^rAd'^rB=NcAB{rA,rB,t)/{VCc), where NcAB{rA,rB,t) is the number of 
triples with A in the volume element d'^rA near the end of the d dimensional radius - vector 
drA and B in the volume element d'^rs near the end of the d dimensional radius - vector 
drB- The analog of Eq. (||) reads 

= {DaVI + DsVl^) Cab {rA, r., t) . (27) 

For the spherically symmetric system, i.e. for Da = Db = D/2 and Ca = Cb, Eq. ( p7[ ) 
reduces to 

dCAB{r,t) ^ j^^i-2d^ f2d-i dCAB{r, t) \ 
dt dr \ dr J 

with the following boundary conditions 

limCAB{r,t)r^^ = CACB, (29) 



DabS; 



2d 



dCAB{r,t) 



dr 



ktCAB{R,t), 



(30) 



r=R 



where r = yr\ + r'^. The solution of the Eq. ( ^Sf ) with the boundary conditions (^91) , ( |30| 
leads to the Eq. (H) for > 1. For = 1. 

Dab 



kth oc 



In 



i?2 • 



(31) 



DABt 



2. Large A and B concentrations (large C — sites) 

Choosing the origin in the center of the C particle, we write the following decoupled 
equations for A and B concentrations 

dCA{B) {r, t) 



dt 

with the boundary conditions 

dCA{B) {r, t) 



DA{B)V^CA(B){r,t) 



(32) 



DA{B)Sd- 



dr 



hCA{r,t)CB{r,t). 



(33) 



When ktCs « AttRcaDa, the solution of the Eqs. dS^^Sp leads to the Eq. (0) for the 
reaction rate with 



kn 



ksmoikt 



ksmol + kfCA 

When ktCB » ^t^RcaDa and ktCA » '^'^RcbDb, it leads to unusual result 



(34) 



dCA 

~dt 



AnCc mm{CARcADA, CbRcbDb)- 



(35) 



IV. DISCUSSION 



Thus, the hoping and diffusion models lead to the similar results for apparent rate con- 
stants — for bimolecular reaction: Eqs. ([TT|) and ( and for threemolecular reactions: 
Eqs. (H, H ) and Eqs. (|l|) — for small C - sites; and Eqs. (|2|, ||)) and Eq. (H, H) — 
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for large C - sites. The diffusion model provides more precise results for the marginal di- 
mensions: d = 2 for bimolecular reactions, Eq. (pA]), and d = 1 for threemolecular reactions, 
Eq. ([3T|). On the other hand, the hoping model provides an easier way for the analysis of 
the correlations structure. Section [111 A| , and for small diffusion coefficient limit, Eq. (pSf). 

In small C — limit the threemolecular ABC correlations are important, since the reaction 
act — annihilation of a triple — changes the concentrations of both A and B in the reaction 
zone from it's maximum value to zero, i.e. S{C^) = 5(C^) = C'^(b) • At the same time, 
since the concentrations are small, the linear boundary condition, Eq. (|30D , is valid. The 
reaction occurs, when A particle joins the CB pair or when B particle joins the CA pair. 
Therefore, the limiting reaction rate, Eqs. (p^ , |3T|) , in this case equals to the rate of jumps 
of the particles A(B) to the reaction zone of CB{A) pair, multiplied by the concentration 
of CB(A) pairs, which meens that kth oc Qksmoi- When Db ^ 0, the reaction is limited 
by the transport of B particles to C and the reaction rate, Eq. (]23|), is the same as for the 
bimolecular reaction B + C ^ with fc^pp =^ ksmoi- 

For large C — sites fluctuations of A and B concentrations in the reaction zone are much 
smaller than the mean values and, therefore, the decoupled (but non linear) equations for A 
and B concentrations, Eqs. ( El) and (B^l, BBf), are valid. For kt — > 0, Eqs. (ER B^l) predict non 



renormalized value, kapp ~ kt\ for kt oo, Eqs. , predict the same reaction 

rate as for diffusion- controlled bimolecular reaction. Note that the reaction "chooses" from 
A and B the reagent with the minimum effective reaction rate.0 More detailed and rigorous 
analysis of the correlator structure will be published in our subsequent publication. 
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